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[i]  A  20-layer,  1/25°  nested  Gulf  of  Mexico  (GoM)  Hybrid  Coordinate  Ocean  Model 
(HYCOM)  has  been  employed  to  examine  the  evolving  three-dimensional  ocean  response 
to  Hurricane  Ivan  during  14-16  September  2004.  Results  from  several  combinations 
of  numerical  experiments  with  and  without  assimilation  of  satellite-altimetry  sea-surface 
height  (SSH)  are  being  analyzed  and  compared  for  the  September  2004  hurricane 
period.  A  comparison  of  simulated  zonal  and  meridional  velocities  using  data  assimilation 
shows  improved  agreement  with  profiler  observations.  The  amplitude  of  the  cold 
wake  (^6°C)  produced  by  these  simulations  compared  reasonably  well  with  the  observed 
changes  in  SST  before  and  after  the  storm;  however,  the  region  of  extreme  cooling 
varied  depending  on  the  simulated  location  of  the  warm  core  eddy  (WCE)  that  had 
detached  from  the  Loop  Current  (LC).  While  the  simulated  location  of  the  WCE  and  LC 
in  the  assimilation  runs  agree  better  with  satellite  altimetry,  the  storm-induced  SST 
cooling  was  40% -50%  greater  than  the  observed  cooling.  Overall,  —64%  of  the  cooling 
was  due  to  vertical  mixing  caused  by  turbulence  generated  from  strong  shear-stress  across 
the  base  of  the  mixed  layer.  Vertical  adveetion  (upwelling)  caused  a  significant  portion 
of  cooling  (23,4%)  in  those  runs  that  included  data  assimilation;  a  three  fold  increase  from 
the  nonassimilative  runs  (7%),  This  enhanced  upper-ocean  cooling  was  caused  primarily 
by  the  p  re  storm  thermal  stratification;  a  shallower  thcrmocline  (—40  m)  and  a  stronger 
upper-thermocline  temperature  gradient  compared  with  the  nonassimilative  rims.  In  all  the 
experiments  the  air-sea  exchange  was  a  small  component  of  the  mixed-layer  heat  budget 
which  overall  accounted  for  —4%. 

Citation:  Prasad,  T.  O.,  and  P.  J.  Hagan  (2007),  Upper-ocean  response  to  Hurricane  Ivan  in  a  1/25“  nested  Gulf  of  Mexico  HYCOM, 
,/.  Geophys.  Res .,  112,  C04013,  doi:10.1029/2006JC0G3695. 


1.  Introduction 

[2]  Warm  water  from  the  Caribbean  Sea  enters  the  Gulf 
of  Mexico  through  the  Yucatan  Straits  and  then  forces  the 
Loup  Current  (LC).  This  anti  eye  Ionic  rotating  current  with 
maximum  Mows  of  1  -2  m  s  1  penetrates  northward  into  the 
Gulf  of  Mexico  and  transports  subtropical  water  with 
markedly  different  temperatures  and  salinities  than  the, 
background  Gulf  of  Mexico  water  between  ocean  temper- 
mures  of  l8°-26CkC.  As  the  LC  intrudes  farther  north ,  the 
instability  processes  associated  with  the  LC  sheds  warm 
core  eddies  (WCE)  having  horizontal  length  scales  of  O 
(200  km  [Elliott,  1982])  at  11-14  month  intervals,  which 

with  speeds  between  and  1  and  14  km 
Crissman,  1986]  over  a  9-12  month 
period  and  eventually  dissipate  along  the  shelf  break  of 
Texas  and  Mexico. 

[3]  Hurricanes  draw  their  energy  from  warm  ocean 
waters.  The  energy  input  from  the  ocean  to  the  hurricanes 
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propagate  westward 
d-1  [Vukovich  and 


and  associated  intensification  depends  largely  on  the  upper- 
ocean  heat  content.  To  quantify  this,  Whitaker  [1967]  first 
used  the  quantity  ‘'hurricane  heat  potential”  which  is 
defined  as  the  vertically  integrated  heat  content  down  to 
the  26°C  isotherm  depth  [Leipper  and  Volgenau t  1972; 
Coni  and  Tvinanes ,  2003],  While  all  hurricanes  attain  their 
maximum  intensity  over  warm  ocean  waters,  there  have 
been  instances  of  sudden  changes  in  intensity  when  passing 
over  areas  of  high  heat  content  [Gon:  and  Trinanes,  2003]. 
In  the  Gulf  of  Mexico,  maximum  hurricane  heat  potential  is 
found  in  the  LC  and  WCE  regions  (~30  x  10*  cal  cm  2), 
where  the  depth  of  26°C  isotherm  may  reach  >100  m 
[Leipper  and  Volgenau,  1972].  These  underlying  preexist¬ 
ing  oceanic  mcsoscale  features  have  far  more  importance  in 
the  heat  and  moisture  fluxes  feeding  die  storm  than  just  SST 
as  noted  in  previous  studies.  Understanding  die  role  of  these 
mcsoscale  features  in  the  intensification  of  hurricanes  is  an 
ongoing  research  topic.  Preliminary  results  have  shown 
their  importance  in  the  sudden  intensification  of  hurricanes, 
Shay  et  al  [2000]  noted  an  abrupt  change  in  the  intensity  of 
hurricane  Opal  (28  September  to  5  October  95)  when  it 
passed  over  a  large  WCE  and  Hong  et  al  [2000]  studied  the 
coupling  between  the  WCE  and  Opal  in  an  atmospheric- 
ocean  coupled  model.  Hurricane  Mitch  occurred  in  the 
Caribbean  Sea  in  October  of  1998  and  Bret  in  the  Gulf  of 
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Mexico  in  August  of  1999  [Coni  and  Tr  in  ones,  2003]  and 
Katrina  in  August  of  2005  [Scharroo  et  ai.,  2005]  also 
experienced  very  similar  intensification  while  passing  over 
LC  and  WCE  regions.  Two  other  hurricanes  that  may  have 
undergone  similar  intensification  are  Ivan  in  September 
2004  and  Rita  in  September  2005,  These  examples  clearly 
demonstrate  the  role  of  warm  oceanic  features  in  providing 
a  positive  feedback  to  the  overlying  storm  by  intensifying 
the  storm, 

[4]  Hurricanes  also  cause  significant  SST  cooling  which 
provides  negative  feedback  to  the  overlying  storm  by 
weakening  the  intensity  of  the  storm.  SST  dropped  as  much 
as  6QC  in  response  to  the  passage  of  Hurricane  Hilda  in 
1964  [Leipper,  1967],  The  lowering  of  SST  during  the 
passage  of  a  hurricane  depends  on  many  factors  such  as 
mixed  layer  depth  (vertical  mixing)  and  thermocline  depth 
(up welling  driven  by  the  wind  stress  curl  of  the  hurricane), 
exchange  of  air-sea  heat  fluxes  and  the  storm's  intensity  and 
translation  speed.  Intense,  slowly  moving  hurricanes  cause 
the  largest  SST  response  [Price,  1981],  When  a  hurricane 
encounters  the  WCE,  die  SST  cooling  may  be  greatly 
suppressed  primarily  due  to  deep  mixed- layer  and  thermo- 
Cline  depths  of  the  WCE. 

[5]  Early  studies  of  upper-ocean  response  to  hurricanes 
included  field  observations  [e*g.,  Leipper,  1967;  Brooks , 
1983;  Sanford  et  ai,  1987;  Shay  and  Etsberry,  1987;  Shay 
et  ai ,  1989,  1998;  Dickey  et  ai ,  1998;  Jacob  et  al. ,  2000] 
and  three-dimensional  numerical  ocean  models  [e*g.,  Price, 
1981;  Price  et  ai,  1994],  The  ocean's  response  to  hurri¬ 
canes  can  be  divided  into  two  stages;  a  forced  and  a 
relaxation  stage  [Price  et  at.,  1994].  In  the  forced  stage, 
hurricane -force  winds  drive  the  mixed-layer  currents,  SST 
cooling  by  vertical  mixing  (entrainment)  and  air-sea  heat 
exchanges  (mainly  due  to  loss  of  latent  heat  flux).  The 
barotropic  response  consists  of  a  geostrophic  current  and  an 
associated  trough  in  sea  surface  height  The  relaxation  stage 
response  following  a  hurricane's  passage  is  primarily  due  to 
inertial-gravity  oscillations  excited  by  the  storm*  The 
mixed-layer  velocity  oscillates  with  a  near-inertial  period 
and  hence  so  does  the  divergence  and  the  associated 
upwdling  and  downwelling. 

[6]  Previous  studies  emphasized  entrainment  mixing  as 
the  dominant  term  in  the  mixed-layer  heat  budget  [Jacob 
ami  Shav,  2003  ].  Jacob  et  ai  [2000]  suggested  that  entrain¬ 
ment  mixing  at  the  base  of  the  mixed  layer  generally 
accounts  for  75% -90%  of  the  cooling  based  on  observa¬ 
tions  while  Price's  [1981]  model  results  indicated  that  85% 
of  die  irreversible  heat  flux  into  the  mixed  layer  was 
through  entrainment  mixing.  Only  about  10%-I5%  of  the 
cooling  in  the  upper  ocean  is  due  to  surface  heat  fluxes  [Price, 
1981]  which  would  range  between  2000  and  3000  W  in-2. 
Jacob  and  Shay's  [2003]  estimate  ranged  from  1 0%  to  30% 
in  the  directly  forced  region,  with  larger  values  in  the  WCE 
and  LC  regions.  Horizontal  advection  is  also  found  to  be 
important  in  the  mixed  layer  heat  balance  during  and  sub¬ 
sequent  to  die  passage  of  hurricanes  [Price,  1981;  Jacob  et 
at.,  2000].  This  contribution  is  particularly  significant  in  the 
eddy  region,  where  maximum  cooling  due  to  geostrophic 
advection  (“0*69°C  d~ 1 )  was  as  large  as  the  surface  heat  flux 
term  in  the  overall  heat  budget  [Jacob  et  ai,  2000].  Also 
the  strong  currents  associated  with  the  preexisting  oceanic 
mesoscale  features  (such  as  warm  and  cold  core  eddies) 


modulate  die  three  dimensional  mixed  layer  heat  budget 
affecting  sensible  and  latent  heat  fluxes  to  the  atmosphere 
through  horizontal  advection, 

[7]  Hurricane  Ivan  (2  24  September  2004)  first  appeared 
on  2  September  2004  as  a  tropical  depression  and  trans¬ 
formed  into  a  category  3  hurricane  after  nearly  three  days. 
Ivan  continued  to  strengthen  as  it  passed  over  the  warm 
waters  of  the  Caribbean  Sea  and  became  a  category  5  hur¬ 
ricane  on  12  September  2004.  During  this  period  the 
atmospheric  pressure  dropped  to  a  minimum  value  of 
910  mb.  Ivan  entered  the  Gulf  on  14  September  2004  and 
began  to  weaken  after  it  passed  the  Loop  Current  region 
and  by  16  September,  had  moved  inland  as  a  category 
3  hurricane.  The  average  translation  speed  of  Ivan  was 
6  m  s-t.  Waves  as  high  as  27.7  m  (91  feet)  were  observed 
when  Hurricane  Ivan  passed  directly  over  the  outer  continental 
shelf  in  the  northeastern  Gulf  of  Mexico  [Wang  et  ai,  2005], 

[ft]  The  oceanic  responses  to  hurricanes  differ  from  one  to 
the  other  in  several  respects  making  the  study  of  the 
processes  difficult.  The  oceanic  response  is  further  compli¬ 
cated  by  preexisting  oceanic  features  such  as  LC  and  WCE 
that  modulate  the  upper  ocean  heat,  mass  and  momentum 
balance  due  to  horizontal  advection  [Jacob  et  al.,  2000; 
Jacob  and  Shay,  2003]*  A  model  that  resolves  these 
mesoscale  features  must  be  used  to  study  the  oceanic 
response  to  hurricanes  in  the  GoM,  Motivated  by  the 
sensitivity  of  the  timing  of  the  eddy-shedding  (WCE) 
behavior  in  the  model  to  changes  in  initial  conditions  which 
arc  vital  for  the  realistic  simulation  of  the  upper-ocean 
response  to  Hurricane  Ivan,  we  carried  out  several  model 
experiments  with  and  without  data  assimilation.  The  model 
is  a  20-layer  1/25°  (^4  km)  nested  Gulf  of  Mexico  Hybrid 
Coordinate  Ocean  Model  (HYCOM).  Two  primary  objec¬ 
tives  of  this  paper  are  (1)  to  assess  the  model's  ability  to 
reproduce  the  observed  behavior  of  the  oceanic  responses 
to  a  hurricane  with  mid  without  data  assimilation  and  (2)  to 
quantify  the  physical  processes  controlling  the  upper-ocean 
thermal  structure  during  Ivan’s  passage  by  evaluating  var¬ 
ious  heat-balance  terms  at  20  m,  50  m  and  100  m  depths. 
The  paper  is  organized  as  follows:  In  section  2,  details  of 
model  configurations,  initial  and  boundary  conditions 
(section  2.1),  surface  forcing  (section  2*2),  vertical  mixing 
(section  2,3)  and  model  experiments  (section  2.4}  are 
presented.  Simulated  quantities  are  compared  with  observa¬ 
tions  in  section  3  followed  by  a  qualitative  description  of  the 
oceanic  response  to  Hurricane  Ivan*  We  quantify  the  upper- 
ocean  physical  processes  governing  the  mixed -layer  heat 
budget  in  section  4*  Results  are  discussed  in  section  5  and 
summarized  in  section  6. 

2,  Model  Description 

[■}]  The  nested-grid  modeling  system  has  a  fine-resolution 
inner  model  embedded  inside  a  coarser-resolution  outer 
model.  The  fine-resolution  nested  GoM  regional  model 
extends  northward  from  18.1°N  and  westward  from 
77. 4°  W.  It  has  a  horizontal  resolution  of  1/25°  (^4  km)p 
twice  the  outer  1/12°  (—8  km)  North  Atlantic  model 
resolution  and  thus  capable  of  resolving  eddies  more 
realistically  and  it  enables  an  even  higher-resolution  nested 
coastal  model  to  be  run.  There  arc  20  hybrid  layers  in  the 
vertical  which  arc  identical  to  the  top  20  layers  in  the  outer 
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2.  M/S 


92°W  8S°W  84*W  9  2*W  88°W  S4°W 


Figure  L  Comparison  of  model  initial  (1  January  2004)  SSH  field  (m)  from  North  Atlantic  1/12°  run 
(left)  with  assimilation  (AS,  ATL-A)  (right)  nonassimilation  (NAS,  ATL-F),  Same  absolute  geostrophic 
velocity  field  from  satellite  altimetry  superimposed  on  both  plots  (cm  s_l)  provides  a  comparison  to 
observed  conditions.  Locations  of  Loop  Current  (LC)  and  WCE  are  different  in  these  runs. 


model  It  should  be  noted  that  the  bottom  5  layers  of  the 
26  layer  outer  model  is  discarded  because  the  densities  of 
the  deepest  5  layers  in  the  outer  model  do  not  exist  in  the 
GoM, 

2.1.  initial  and  Boundary  Conditions 

[m]  The  lateral  boundary  conditions  for  the  GoM  nested 
model  come  from  the  1/1 2°  (^8  km)  North  Atlantic  data- 
assimilative  HYCOM  system  (ATL-A),  In  the  nested  model 
most  of  the  How  enters  the  domain  on  the  southern  and 
southeastern  boundaries  and  exits  through  ihe  Florida  Strait. 
At  the  open  boundaries,  buffer  (or  boundary  relaxation) 
zones  typically  less  than  one  degree  wide  (where  the  two 
model  grids  overlap)  are  used  to  relax  the  baroclinic  model 
temperature,  salinity,  pressure  and  velocity  components 
once  per  day  towards  the  outer  model  fields  (ATL-A)  with 
an  e-folding  time  scale  of  1  to  10  days.  This  procedure  has 
proven  to  be  very  robust  The  method  of  characteristics 
[Browning  ami  Kreiss,  1982,  1986]  is  used  for  die  baro- 
tropic  open  boundary  conditions  on  velocity  and  pressure. 
HYCOM  currently  has  a  robust  capability  for  nesting  with 
other  HYCOM  grids  with  similar  vertical  design.  With  this 
method  all  information  passes  from  the  coarse  outer  grid  to 
the  finer  inner  grid  and  does  so  offline,  meaning  that  the 
nested  model  does  not  run  concurrently  with  the  outer 
model. 

[n]  The  coarse -resolution  outer  model  (ATL-A)  extends 
from  28°S  to  7G°N  including  the  Mediterranean  Sea  (see 
Chassignet  et  at.  [2007]  for  details).  The  vertical  resolution 
consists  of  26  hybrid  layers,  with  the  top  layer  typically  at 
its  minimum  thickness  of  3  nr  The  bathymetry  used  in  the 
model  is  derived  from  World  Ocean  Elevation  data 
(ETOP05)  topography.  The  northern  and  southern  bound¬ 
aries  are  treated  as  closed,  but  are  outfitted  with  3°  buffer 
zones  in  which  temperature,  salinity  and  pressure  are 
linearly  relaxed  toward  their  seasonally  varying  climatolog¬ 
ical  values  thereby  simulating  the  thermoh aline  overturning 
circulation.  The  data  assimilation  in  the  1/1 2°  North  Atlantic 
HYCOM  system  consists  of  assimilating  daily  operational 
Modular  Ocean  Data  Assimilation  System  (MODAS)  1/4° 
SSII  analysis  of  available  real-time  satellite  altimeter  obser¬ 
vations,  The  Cooper  and  Haines  [1996]  technique  is  used  to 


project  the  surface  information  to  the  interior  of  the  ocean. 
Relaxation  to  the  MODAS  SST  analysis  derived  from  the 
5-channel  Advanced  Very  High  Resolution  Radiometers 
(AVHRR)  is  also  included.  A  similar  configuration  of  tile 
North  Atlantic  HYCOM  without  data  assimilation  (ATL-F) 
has  been  running  in  parallel  with  the  assimilative  system. 

[12]  The  initial  conditions  for  the  simulations  performed 
here  are  extracted  from  these  two  configurations  of  the 
North  Atlantic  HYCOM  (ATL-A  and  ATL-F).  All  mode! 
integrations  are  started  on  January  1,  2004,  Figure  1  shows 
the  initial  SSH  (1  January  2004)  from  assimilative  (ATL-A) 
and  nonassimilative  (ATL-F)  Atlantic  1/12°  model  runs.  For 
comparison  to  observed  conditions,  the  same  geostrophic 
velocity  field  from  merged  satellite  altimetry  products  is 
superimposed  on  both  plots.  As  expected,  the  location  of 
WCE  and  the  northward  penetration  of  LC  in  these  two 
model  configurations  are  different  as  indicated  by  the  SSH, 
The  location  of  the  WCE  and  LC  in  the  assimilation  run 
agrees  witli  the  observations  while  the  northern  edge  of  LC 
in  the  free  run  occurs  farther  south  and  there  is  no  indication 
of  WCE.  It  is  these  differences  that  subsequently  lead  to 
differences  in  the  WCE  location  and  SST  cooling  during 
Ivan’s  passage, 

2.2.  Surface  Forcing 

[13]  The  model  is  driven  by  fields  of  10  m  wind  speed, 
vector  wind  stress,  2  m  air  temperature,  2  m  atmospheric 
humidity,  surface  shortwave  and  long -wave  heat  fluxes,  and 
precipitation.  These  fields  are  extracted  from  three-hourly 
1°  horizontal  resolution  Navy  Operational  Global  Atmo¬ 
spheric  Prediction  System  (NOGAPS)  reanalysis  product 
for  the  year  2004.  Surface  latent  and  sensible  heat  tluxes, 
along  with  evaporation  arc  calculated  using  bulk  formulae 
during  the  model  run  time  using  model  SST.  Tins  has  an 
implied  restoring  term,  pulling  the  model  produced  SST 
towards  the  specified  air-temperature  thereby  minimizing 
model  SST  drift. 

[14]  The  NOGAPS  wind  stress  is  computed  using  the 
drag  formulation  from  Louis  [1979].  HYCOM  includes 
several  air-sea  flux  parameterization  schemes  to  determine 
Ihe  exchange  coefficients  for  heat  fluxes.  While  the  values 
of  these  coefficients  under  high  winds  are  an  ongoing 
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research  topic  [Powell  et  al. ,  2003],  the  exchange  coeffi¬ 
cients  for  heat  (i.e.,  CL ,  and  Cs)  used  here  have  a  simple 
polynomial  dependence  on  wind  speed,  and  a  linear  depen¬ 
dence  on  the  air-sea  temperature  difference  {see  Kara  et  al 
[2002]  for  details).  Briefly,  latent  heat  flux  { Q /  }  and  sensible 
heat  flux  (Qs)  are  computed  using  the  equations 

Qs^CsCpMZ'-Ts), 

Ql  —  ^ i L pa ( <7 u  —  lh)> 

CL  =  Qq  +  Cu(Ts  -  Ta),  Cs  =  0.95Qt 

Q0  =  1  O'3  [0.8 1 95  +  0,0506  Va  -  0.0009  (K,)2J, 

Ct|  =  10“,l-0.0154  + 0.5698(1/^)  -0.6743(1/^,)^, 

Vu  is  wind  speed  at  10  m  above  sea  level;  Ta  is  air- 
temperature;  7V  is  model  SST;  Cp  is  specific  heat  of  air 
(1004.5  3  kg  1  K  ');  and  L  is  latent  heat  of  vaporization 
{2,5  x  10fl  J  kg-1).  The  air-density  (kg  m-3)  is  determined 
using  the  ideal  gas  law;  p„  =  100  PJ[R^„y(Ta  +  273.16)j, 
where  Pa  is  set  to  L  0  L  3  liPa.  The  mixing  ratio  values  for  air 
(tja  at  Tit)  and  sea  (qs  at  7])  are  calculated  using  a  simplified 
version  of  the  original  formulation  for  saturation  vapour 
pressure  (es)  presented  by  Buck  [1981], 

23,  Vertical  Mixing 

[\5]  HYCOM  includes  several  turbulence  closures  to 
represent  vcrtieal  mixing.  In  addition  to  the  Kraus-Tumer 
model  used  in  MICOM  and  modified  to  run  wiih  hybrid 
coordinates,  IIYCOM  has  been  equipped  with  the  Price- 
Wcllcr-Pinkel  dynamical  instability  model  (PWP  [Price 
et  al,  1986]),  the  K- Profile  Parameterization  (KPP  [Large 
et  alf  1994,  1997]),  the  Mdlor-Yamada  level  2,5  (MY 
[Me  I  hr  and  Yam  ad  a,  1982])  and  the  NASA  Goddard 
Institute  for  Space  Studies  level  2  (GISS  [Canuto  et  al , 
2001,  2002]).  An  evaluation  of  these  vertical  mixing 
schemes  in  the  IIYCOM  in  low -resolution  climatological 
simulations  of  the  Atlantic  Ocean  is  discussed  by  HaUiwell 
[2004].  Technical  details  of  the  implementation  of  these 
vertical  mixing  schemes  have  been  included  in  the  HYCOM 
User's  Manual  [Blech  et  al. ,  2002].  The  experiments 
reported  here  included  the  G1SS  mixing  scheme.  This 
vertical  mixing  model  was  constructed  using  the  Reynolds 
stress  model  (sec  Canuto  et  al  [2001,  2002,  2004]  for 
details).  The  GISS  model  included  a  more  physical  repre¬ 
sentation  of  several  aspects  of  the  Reynolds  stress  model,  in 
particular  a  better  representation  of  the  velocity  and  tem¬ 
perature  pressure  correlations. 

[ig]  The  simulation  of  mixed  layer  deepening  and  cooling 
during  the  storm's  passage  strongly  depends  on  the  choice 
of  a  vertical  mixing  scheme.  Recent  studies  of  uppcr-oecan 
responses  to  Hurricane  Gilbert  based  on  several  bulk  mixed 
layer  entrainment  schemes  revealed  significant  differences 
in  the  heat  and  mass  budgets  [Jacob  et  al,  2000;  Jacob  and 
Shay ,  2003].  Though  most  schemes  differ  in  their  physical 
basis,  direct  comparisons  indicate  no  clear  advantage  bet¬ 
ween  several  schemes.  While  these  discrepancies  among  the 
schemes  make  it  hard  to  justify  the  suitability  of  a  particular 
scheme,  GISS  has  a  number  of  advantages  over  other 


schemes  including  (1)  the  higher  predictability  of  die 
turbulent  mixing  Rt  ^  0(1)  compared  to  MY-type  models 
which  cut  off  mixing  too  early  leading  to  shallow  mixed 
layers  [Martin,  1985]  (2)  the  number  of  adjustable  para¬ 
meters  is  one  {compared  to  11  in  KPP)  and  (3)  the 
latitudinal  dependence  of  the  diffusivities  in  the  thcrmocline 
(which  KPP  does  not  account  for).  Thus  KPP  is  basically  a 
diagnostic  model  whereas  GISS  is  prognostic  and  thus 
much  more  suitable  for  climate  studies.  The  primary  con¬ 
clusions  drawn  here  using  GISS  also  stand  valid  for  other 
vertical  mixing  schemes. 

2,4,  Model  Experiments 

[it]  The  major  mesoscale  features  of  oceanic  circulation 
in  the  GoM  arc  ihc  LC  and  WCE  that  detach  Ifom  it. 
During  its  passage.  Hurricane  Ivan  encountered  both  of 
these  predominantly  ant  icy  clonic  circulation  features. 
Therefore,  realistic  prestorm  conditions  including  the  major 
currents  and  eddies  prior  to  the  passage  of  Ivan,  arc  vital  for 
the  accurate  simulation  of  the  upper-ocean  response.  This  is 
achieved  here  by  performing  a  series  of  model  experiments 
using  a  nested  GoM  regional  model.  For  simplicity,  our 
control  (main)  run  does  not  include  data  assimilation.  The 
initial  condition  for  this  control  run  comes  from  a  non- 
assimilarive  version  of  1/12°  North  Atlantic  system 
(ATL-F)  and  the  simulation  does  not  include  data  assimi¬ 
lation  (free  run).  As  we  shall  see  in  the  following  sections, 
the  control  run  failed  to  simulate  the  WCE  location  con¬ 
sistent  with  the  observations  during  the  Hurricane  Ivan 
period.  To  resolve  this  issue,  we  design  four  additional 
experiments  by  (1)  changing  the  initial  conditions  and  (2) 
including  the  data  assimilation  during  the  model  run.  In  the 
first  experiment,  we  initialize  the  control  run  with  more 
realistic  initial  conditions  from  an  assimilation  version  of 
1/12°  North  Atlantic  model  (ATL-A).  Hie  second  experi¬ 
ment  is  an  exact  repeat  of  the  control  run  except  that  the 
simulation  included  data  assimilation  throughout  the  period. 
Comparison  of  this  experiment  with  the  control  run  will 
determine  the  impact  of  the  assimilation.  We  devise  a 
scheme  to  identify  these  experiments  by  taking  into  account 
two  aspects  of  the  simulation  (1)  initial  conditions  and 
(2)  run-lime  specifications.  On  the  basis  of  the  model's 
initial  condition  (e.g.  assimilation,  free-run)  and  run-time 
specifications  (e.g.  assimilation,  free-run),  we  denote  these 
simulations  FF  (control  run),  AF  (experiment  1>  and 
FA  (experiment  2)  respectively  (e.g.  in  AF,  A  stands  for 
assimilation,  the  initial  condition,  and  F  stands  for  tree-run, 
also  nonassimilation,  the  run-time  specification).  An  exact 
repeat  of  AF  but  with  data  assimilation  (AA)  produced 
identical  results  with  that  of  FA  and  therefore  is  not 
presented  here. 

[is]  Two  additional  experiments  that  arc  variants  of  these 
three  runs  are  also  performed  fo  address  specific  questions 
on  the  effectiveness  of  the  non  assimilative  versus  assimi¬ 
lative  runs.  The  first  experiment  is  an  exact  repeat  of  the 
control  run  except  that  live  assimilation  is  applied  only 
through  1  —30  September  2004  (FF-A9).  FF-A9  differs 
from  FA  in  that  assimilation  is  applied  from  I -January  to 
31 -December  2004.  Ln  the  second  experiment,  the  control 
run  is  repeated  with  data  assimilation  except  during  the 
period  1  -30  September  2004  (FA-F9).  These  two  simula¬ 
tions  are  carried  out  only  for  live  period  1-30  September 
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Table  1.  Model  Experiments 


Experiments 

LC." 

LC. 

Time 

1C. 

Experiments 

Run 

Period 

FF 

NASh 

1-Jan  04 

ATL-F 

NAS 

Jan-Dcc 

AF 

AS" 

LJan  04 

ATL-A 

NAS 

Jon-Dee 

FA 

NAS 

l-Jan  04 

ATL-F 

AS 

Jan-Dec 

FF-A9 

NAS 

1  -Sep  04 

FF 

AS 

1-30  Sep 

FA-F9 

AS 

I -Sep  04 

FA 

NAS 

1-30  Sep 

initial  conditions. 
^NAS,  non  assimilation. 
CAS,  ussimilation. 


2004  by  taking  initial  conditions  (J -September  2004)  from 
the  experiments  FF  and  FA,  respectively.  Thus,  in  FF-A9, 
FF  denotes  the  baseline  run  from  which  initial  conditions 
are  taken  and  A9  indicates  the  run-time  specifications 
(assimilation  during  1-30  September  2004).  The  initial 
conditions  (l -September  2004)  for  FA-F9,  like  AF,  includes 
the  correct  location  and  structure  of  the  VVCE  and  LC. 
These  model  experiments  arc  summarized  in  Table  l.  It 
should  be  noted  that  in  all  the  experiments  reported  here, 
SST  evolves  as  a  result  of  physical  processes  without 
relaxing  il  to  MODAS  SST  This  enables  us  to  compare 
the  mixed  layer  heat-budget  with  and  without  data  assim¬ 
ilation.  The  data  from  these  experiments  analyzed  during 
the  Ivan  period  (September,  2004)  are  presented  in  the 
following  sections. 


3*  Results 

[19]  Results  from  three  model  experiments  are  presented. 
The  impact  of  data  assimilation  on  the  model  simulations 
are  assessed  by  comparing  SSH,  SST  and  velocity  fields 
from  FF  and  FA  with  observations.  The  SSH  variance  (nr) 
from  FF  and  FA  is  compared  with  that  derived  from 
altimetry  for  the  year  2004  in  Figure  2.  Daily  SSH  is  used, 
Tire  energy  level  of  the  WCE  is  realistic  and  the  eddy 
detachment  from  the  LC  occurs  at  approximately  the 
correct  latitude.  The  location  of  extreme  variance  coincides 
with  the  WCE  and  its  westward  propagation.  The  amplitude 
is  larger  (smaller)  in  the  simulation  without  (with)  data 
assimilation. 

[zq]  Time  series  of  SST  from  National  Data  Buoy  Centre 
(NDBC)  buoy  42001  located  on  the  left  side  of  Ivan's 
track  (89.66° W,  25.84*N)  is  compared  with  model  SST 
{Figure  3a).  The  inset  plot  is  for  the  September  time-period; 
the  focus  of  the  study.  It  should  be  noted  that  the  model 
SSTs  arc  the  mixed  layer  temperature,  while  the  buoy 
measures  the  actual  surface  temperature.  Simulated  SSTs 
between  May  and  August  2004  are  in  good  agreement  with 
ihe  buoy  SSTs  while  they  depart  systematically  during  the 
poststorm  period.  In  September  the  buoy  SST  indicated  two 
periods  of  cooling;  the  first  event  took  place  during  the 
storm  and  the  second  major  cooling  occurred  after  Ivan 
made  landfall.  The  cooling  during  the  storm  was  primarily 
due  to  the  net  surface  heat  loss  from  the  ocean  (—300  W  m_:f) 
and  wind -driven  turbulent  mixing  due  to  strong  wind 
(16  m  s  ,  figures  not  shown).  The  cooling  during  the 
poststorm  period  was  not  driven  by  surface  heat  flux,  but 
was  mosi  likely  due  to  adveetion  or  upwelling. 

[21]  Mixed-layer  zonal  and  meridional  velocity  compo¬ 
nents  from  FF  and  FA  are  compared  with  the  altimetry 


derived  absolute  geostrophic  velocity  vectors.  Mindful 
that  the  model  velocity  fields  include  a  wind- driven  com¬ 
ponent,  the  one-to-one  comparison  should  be  interpreted 
carefully.  The  area  averaged  (88°  86° W,  24°  26°N)  ve¬ 
locity  fields  are  shown  in  Figures  3b  and  3c.  The  velocity 
components  from  FA  agree  favorably  with  the  altimetry, 
but  they  differ  significantly  with  FF.  The  large  differences 
in  the  velocity  components  that  occurred  during  April- 
August  are  associated  with  the  differences  in  the  location  of 
the  WCE  among  the  two  model  runs  which  is  reflected  in 
SSH  (Figure  3d).  The  separation  of  WCE  from  the  LC  in 
FF  took  place  in  April  (indicated  by  peak  SSH)  whereas 
that  in  FA  occurred  in  August  consistent  with  observations. 
As  we  shall  discuss  in  the  following  sections,  the  different 
time  of  separation  of  WCE  from  the  LC  in  these  runs  led  to 
changes  in  ihe  advective  heal  transport  during  the  passage 
of  Ivan. 

[22]  Ivan  passed  directly  over  an  array  of  14  Acoustic 
Doppler  Current  Profilers  (ADCPs)  deployed  along  the 
outer  continental  shelf  and  upper  slope  in  the  northeastern 
Gulf  of  Mexico  as  part  of  the  Slope  to  Shelf  Energetics  and 
Exchange  Dynamics  (SEED)  project  [Wang  et  aLt  2005] 
(W,  J.  Teague  et  aL,  Observed  oceanic  response  over  ihe 
upper  continental  slope  and  outer  shelf  during  Hurricane 
Ivan,  submitted  to  Journal  of  Physical  Oceanography* 
2006,  hereinafter  referred  to  as  Teague  et  aL,  submitted 
manuscript,  2006).  Simulated  current  profiles  were  com¬ 
pared  with  ADCP  measured  zonal  (U.  m  s  ')  and  merid¬ 
ional  (K  m  s^1)  ocean  currents  from  moorings  M3  and  M 1 3 
situated  to  the  right  of  Ivan’s  track  to  evaluate  the  models 
skill  in  predicting  the  ocean  response.  The  choice  of  these 
moorings  for  comparison  was  made  according  to  their 
locations;  M3  being  located  in  the  outer  shelf  (87.84°W, 
29.47°N)  at  a  depth  of  ^60  m  and  MI3  being  located 
along  the  continental  slope  (87.83DW,  29.16°N)  at  a  depth 
of  1000  m,  thereby  enabling  us  to  compare  the  ocean's 
response  to  Ivan  across  the  shelf  break.  M3  (Ml  3)  recorded 
current  profiles  with  2  m  (10  m)  vertical  resolution  every  15 
(60)  minutes  with  an  accuracy  of  0,5%  ±  0.5  cm  s~* 
('-l  cm  s~‘).  A  24  hour  Parzen  window  smoothing  has 
been  applied  on  these  data  for  direct  comparison  with 
model  fields.  A  detailed,  nonsmooihed,  version  of  the 
velocity  components  from  all  14  moorings  are  plotted  and 
discussed  by  Teague  et  al.  (submitted  manuscript,  2006), 

[23]  Time  series  of  zonal  and  meridional  velocity  com¬ 
ponents  from  moorings  M3  and  Ml 3  in  the  upper  50  and 
500  in,  respectively,  and  the  corresponding  fields  from 
model  simulations  FF  and  FA  are  depicted  in  Figure  4. 
Overall,  simulated  currents  compared  well  with  the  obser¬ 
vations  and  zonal  component  dominated  the  fiow  in  part 
due  to  the  flow  being  parallel  to  the  isobaths.  In  the  outer 
shelf  (M3),  the  zonal  flow  during  the  prestorm  period 
(5-12  September)  was  eastward  with  peak  velocities  of 
0,5  ms3  at  20-3(1  m.  Both  simulations  reproduced  this 
flow  with  assimilation  run  FA  being  the  closest  match.  The 
storm-induced  flow  reversal  in  the  upper  50  111  was  evident 
in  both  the  observations  and  simulations  with  maximum 
velocities  exceeding  1  m  s_l.  The  subsurface  intensified 
eastward  fiow  (0.4  m  s”1)  at  M3  seen  after  25  September 
was  less  pronounced  in  the  model  (0.2  m  s_1),  In  the  wake 
of  Ihe  storm,  the  zonal  velocity  oscillated  with  alternating 
eastward  and  westward  motion  indicating  near-inertial 
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currents  which  the  model  successfully  simulated.  Though 
much  weaker  than  V ,  the  simulated  meridional  velocity 
agreed  reasonably  well  with  ADCP  with  one  major  excep¬ 
tion;  a  strong  southward  How  (0,3  m  s-1)  seen  in  ADCP 
data  between  17  and  21  September  was  absent  in  the 
model 

[24]  The  upper-ocean  currents  response  to  Ivan  at  Ml 3 
was  not  significantly  different  from  that  at  M3.  The  storm- 
forced  maximum  simulated  zonal  velocities  in  FF  slightly 
underestimated  (03  m  s  ’}  the  observed  maximum  of 
0.5  m  s" 1  below  —200  m,  but  agreed  better  in  FA.  Prestorm 
and  poststorm  velocities  were  also  greatly  improved  when 
simulation  included  the  data  assimilation.  For  example,  the 
magnitude  of  the  prestorm  eastward  flow  in  FF  was  30%- 
40%  huger  than  the  ADCP  data,  but  agreed  better  with  FA, 
The  agreement  between  the  observed  and  modeled  merid¬ 
ional  velocities  was  reasonably  good.  A  similar  degree  of 
agreement  has  also  been  evident  for  all  1 2  other  mooring 
locations  which  are  not  discussed  here. 

[25]  Time  series  of  zonal  and  meridional  wind  stresses 
from  NOGAPS  are  shown  in  Figure  5a.  A  maximum  wind 
stress  of  1.82  N  m~“  was  attained  during  the  passage  of 
Ivan,  The  time  series  of  mixed  layer  temperature  (SST), 
thickness  (MLD),  salinity  (SSS),  and  surface  heat  flux 
(SHF)  averaged  for  the  region  88°-86°Wt  24°-26°N  from 
these  three  experiments  are  shown  in  Figures  5b-  5e  and 
Tropical  Mierowave/lmugcr  (TMI)  derived  SST  for  com¬ 
parison  in  Figure  5b.  In  general,  the  SST  cooling  during  the 
stomi  was  in  reasonable  agreement  with  the  observations. 
The  prestorm  MLD  of  a  M0  m  increased  to  —45  m  during 
the  storm  and  the  corresponding  SSTs  or  mixed  layer 
temperature  decreased  from  —28 ,5°  to  — 25°C  (— 3.5dC 
decrease)  except  in  the  run  that  included  assimilation.  The 
cooling  simulated  by  FF  was  in  better  agreement  with  the 
TMI.  Despite  a  shallower  ML  (—30  m)t  the  SST  cooling  in 
FA  was  far  more  intense  {— 22.5°C)  than  the  other  two 
experiments  and  TMI  It  is  likely  that  part  of  the  increase  in 
the  MLD  is  offset  by  the  up  welling  which  reduced  the  mixed 
layer  thickness  by  15  m  (discussed  further  in  section  4). 
Greatbatch  [1985]  noted  the  role  played  by  the  upwdling  of 
water  and  the  associated  reduction  in  mixed-layer  depth 


in  lowering  SST  during  the  passage  of  a  storm.  Poststorm 
model  SSTs  were  cooler  by  — 2°-  4°C  than  TMI  These 
SST  differences  among  the  simulations  were  reflected  in  the 
net  surface  heat  flux  through  surface  latent  and  sensible  heat 
fluxes  terms,  which  were  calculated  using  the  model  SST. 
The  net  surface  heal  flux  showed  a  loss  of  —150  W  m-2 
to  the  atmosphere  during  the  storm  which  partly  accounted 
for  SST  cooling  in  A F  and  FF,  but  it  was  only  about 
50  W  m_1  in  FA.  As  expected,  this  large  heat  toss  was 
primarily  due  to  the  high  latent  heat  flux  exchange  caused 
by  the  high  winds,  Poststorm  MLD  in  all  experiments  dra¬ 
matically  dropped  to  —5  m  with  SSTs  steadily  increasing. 
The  poststorm  wanning  was  in  part  a  result  of  weaker  winds 
along  with  cool  SSTs  which  reduced  the  latent  heat  loss 
significantly  leading  to  150  W  m-2  heat  flux  gain  by  the 
ocean.  What  is  important  here  is  that  surface  heat  loss  to  the 
atmosphere  was  not  a  significant  contributor  to  SST  cooling 
during  the  storm. 

[26]  Consistent  with  SST  changes,  SSS  also  showed 
significant  differences  among  the  three  experiments 
though  freshwater  fluxes  (evaporation  -  precipitation) 
from  these  experiments  were  nearly  the  same.  The  max¬ 
imum  surface  freshening  occurred  in  FA  where  it  dropped 
to  35.6  psu  from  prestorm  value  of  35.8  psu.  The  lowest 
SSS  associated  with  the  storm  occurred  one  day  later  than 
the  SST  minimum.  In  the  Gulf  of  Mexico,  low  salinity 
waters  are  located  below  —100  m  (except  in  the  northern 
Gulf)  with  salinity  increasing  with  depth  in  the  upper 
100  m.  If  this  surface  freshening  occurs  as  a  result  of 
upwelling,  then  it  has  to  come  from  much  deeper  depths. 
On  the  other  hand,  upwelling  or  mixing  of  cooler  water 
from  a  relatively  shallow  thermocline  yielded  instant  SST 
cooling. 

[27]  The  complex  nature  of  interactions  between  the 
WCE  and  Hurricane  Ivan  are  evident  in  Figure  6.  TMI 
derived  SST  (shaded)  superimposed  with  altimeter  SSH 
(contours)  before  ( 1 0  September  2004)  and  after  (16  September 
2004)  the  passage  of  Hurricane  Ivan  (SSTs  are  three  day 
averages  ending  on  10  and  1 6  September  2004,  respectively) 
and  the  corresponding  model-derived  SST  (snapshot)  and 
SSH  from  the  three  experiments  are  displayed  in  Figure  6 
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Figure  3-  Simulated  SSTs  from  FF  (without  assimilation)  and  FA  {with  assimilation)  is  compared  with  a 
National  Data  Buoy  Centre  (NDBC)  buoy  (42001)  at  89.66°W,  25.84°N  located  to  the  west  of  Ivan’s 
track.  The  inset  plot  is  for  September  time  frame;  the  focus  of  the  study.  Comparison  of  model  (b)  zonal 
and  (c)  meridional  velocity  components  (ms  ')  with  observed  absolute  geostrophic  velocities  (altimetry) 
and  (d)  SSH  (m).  These  fields  are  averaged  over  the  2°  x  2°  box  24°-26*N)  shown  in 

Figure  6.  A  five-day  smoothing  has  been  applied  to  daily  model  fields  to  make  them  consistent  with  five- 
day  altimetry  data.  The  vertical  lines  in  each  plot  indicate  the  time  of  arrival  and  landfall  of  Ivan  in  the 
Gulf  of  Mexico, 


The  right  panel  of  Figure  6  shows  the  corresponding  SST 
differences  (ASST).  Prestorm  SST  (10  September  2004)  in 
the  Gulf  was  >30°C  almost  everywhere  resulting  from  the 
intense  summer  solar  heating  except  a  small  region  south¬ 
west  of  the  Florida  coast  (remnants  of  cooling  from  Hurri¬ 
cane  Frances  about  two  weeks  earlier).  The  altimeter  SSH 
revealed  the  existence  of  a  cyclonic  cold-core  eddy  (CCE, 
-87°W,  25°N)  and  an  anticyclonic  WCE  (-89.2°W, 
26.7°N)  along  Ivan’s  track,  P  re  storm  SST,  however,  did 
not  indicate  the  presence  of  these  mesoscale  eddies.  What  is 
interesting  here  is  that  the  amplitude  of  SST  cooling  in  the 
wake  of  the  hurricane  was  not  continuous.  The  maximum 
(minimum)  cooling  occurred  outside  (inside)  the  WCE  and 
LC  regions  suggesting  that  the  upper-ocean  cooling  was 
directly  influenced  by  these  preexisting  mesoscale  features. 
In  contrast  to  the  well-documented  rightward  bias  of  the 
SST  cooling  response  to  most  hurricanes,  the  maximum 
SST  cooling  occurred  to  the  left  of  the  Ivan’s  track, 
especially  west  of  the  LC  region  (TMI  SST  during 
September  16),  Leipper's  [1967]  observations  of  the  SST 


response  to  Hurricane  Hilda  did  not  show  a  rightward  bias 
(see  their  Figure  8), 

[zk]  Though  hurricane-induced  SST  cooling  was  certainly 
evident  in  all  three  experiments,  the  details  of  the  cooling 
especially  the  location  of  maximum  cooling  differed  from 
the  TMI  and  among  these  three  simulations.  These  differ¬ 
ences  were  caused  primarily  by  the  differences  in  the 
location  of  warm  and  cold  core  eddies  which  are  clear 
from  the  superimposed  SSH  fields.  While  LC  and  WCE  in 
AF  and  FA  were  located  in  the  vicinity  of  the  observed 
location,  the  LC  in  FF  penetrated  farther  north  and  WCE 
was  situated  ^395.7  km  farther  southwest  (92.2°W, 
26. 1°N)  than  altimetry  (SS.3°W,  26  8°N),  The  CCE  in  FF 
was  located  slightly  (~78  km)  farther  northwest  (87r5°W, 
25,5°N)  than  altimetry  (87. 1°W,  24,9°N).  Consequently,  the 
region  of  maximum  cooling  occurred  farther  north  of  the 
observed  minimum.  Coincident  with  the  southward  shift  of 
the  WCE,  the  maximum  SST  cooling  in  AF  occurred  farther 
south  compared  with  TMI,  The  location  of  the  greatest 
cooling  in  FA  agreed  with  the  TMI  observations,  but  the 
amplitude  of  simulated  SST  differences  (ASST)  induced  by 
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Figure  4*  Time  series  of  observed  zonal  {( J)  and  meridional  (F)  velocity  components  fiom  moorings 
M3  (87.84°W,  29.47°N)  and  M13  (87T83°W,  29.16°N)  in  the  upper  50  and  500  m,  respectively,  and  the 
corresponding  fields  from  model  simulations  FF  and  FA.  A  24 -hour  Parzcn  window  smoothing  has  been 
applied  on  these  mooring  data.  Contour  interval  is  0.1  m  s“\  These  moorings  are  located  to  the  right  of 
Ivan's  tracks  Daily  model  fields  are  used*  Note  different  color  scaling  is  used  for  U  and  F  Positive  U  (F) 
indicates  eastward  (northward)  flow. 


8  of  L8 


C04G13 


PRASAD  AND  HOGAN:  MIXED  LAYER  RESPONSE  TO  HURRICANE  IVAN 


CO-JOB 


(88°-86°W,  24°-26°N) 


3Q.Q 

23.0 

26.0 

m 

24  0  LO 
220 


300 

160 


0  5 


-150 

-300 

6 


X 

</l 


CL 

I 

LU 


0 


Figure  5,  Time  series  of  (a)  zonal  (east-west)  and  meridional  (north-south)  wind  stress  (N  m~2)t  model 
(b)  mixed  layer  temperature  (SST,  °C),  (c)  mixed  layer  depth  (MLD,  m),  (d)  surface  heat  flux  (SHF, 
W  m  2X  and  (e)  mixed  layer  salinity  (S$$,  psu)  from  AF.  FF  and  FA.  Also  included  in  Figure  5b  is  the 
TM1  derived  SST,  and  model  evaporation-precipitation  (E-P,  kg  m"3  s->)  is  included  in  Figure  5e. 


Ivan  was  as  high  as  10°C  and  was  spread  over  a  much 
larger  area. 

[24]  In  contrast  to  TML  the  signature  of  the  LC  and  WCE 
were  clearly  identifiable  in  the  prestorm  SST  maps  in  all 
experiments.  Simulated  prestorm  SSTs  (especially  north¬ 
west  of  LC)  were  generally  colder  than  those  derived  from 
the  TMI  by  M°C.  The  coincidence  of  the  extreme  cooling 
within  the  cyclonic  eddy  in  all  experiments  during  Ivan’s 
passage  suggests  that  upwelling  was  likely  an  important 
mechanism  for  upper-ocean  cooling.  The  poststorm  SSH 
dropped  from  the  prestorm  conditions  by  MO  cm  due  to 
shoaling  of  upper  thermocline  driven  by  the  curl  of  the 
wind  stress.  While  the  amplitude  of  simulated  SST  differ¬ 
ences  were  within  the  observed  limits  (except  for  FA),  there 
were  significant  differences  between  the  simulated  and 
observed  location  of  maximum  cooling.  What  is  however 
interesting  in  both  simulated  and  observed  cooling  is  that 
the  cooling  was  only  about  1°-2°C  within  WCE  and  LC 
regions.  These  results  suggested  a  clear  modulation  of 
upper-ocean  cooling  by  the  WCE  and  LC  systems, 

[30]  The  background  circulation  associated  with  these 
systems  combined  with  directly  forced  wind-driven  cur¬ 
rents  during  the  passage  of  the  storm  also  played  an 


important  role  in  the  modulation  of  SST  cooling  via 
horizontal  adveetion.  To  demonstrate  this,  we  show  pres- 
lorm  (10  September)  and  storm  (15  September)  ocean 
surface  currents  from  AF,  FF  and  FA  in  Figure  7,  Prestorm 
surface  circulation  in  the  Gulf  was  largely  geostrophic;  the 
magnitude  of  the  LC  {—1  m  s_l)  is  comparable  with 
the  absolute  geostrophic  velocity  vectors  derived  from  the 
satellite  altimetry  (Figure  t).  Because  of  die  westward 
displaced  WCE  in  FF,  the  LC  penetrated  farther  north 
than  in  the  other  two  experiments.  During  the  storm, 
directly  forced  wind-driven  currents  dominated  the  surface 
circulation  in  all  experiments  and  the  structure  of  the  LC 
was  not  disoemabic.  However,  the  effect  of  prestorm 
circulation  associated  with  the  WCE  and  LC  on  the  storm 
induced  currents  was  evident.  For  example,  the  westward 
flowing  current  to  the  southern  flank  of  the  WCE  in  AF 
was  slightly  stronger  than  that  in  FF.  The  strongest 
currents  exceeding  2ms"1  were  found  to  the  right  of 
Ivan’s  tTack,  consistent  with  the  most  intense  winds  asso¬ 
ciated  with  the  eye-wall  being  located  in  the  northeast 
quadrant  of  the  hurricane,  leading  to  large  current  shear 
across  the  mixed  layer.  To  the  left  of  the  track  the  currents 
were  somewhat  weaker.  This  suggests  that  the  horizontal 
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Figure  6*  Prestorm  ( 10  September)  and  poststonn  (16  September)  SSTs  and  die  corresponding  prestorm 
and  poststonn  SST  differences,  ASST  (shaded,  °C)  from  the  Tropical  Microwave/Imager  (TMI)  and  from 
model  experiments  AF  FF  and  FA  arranged  from  top  to  bottom  panels,  respectively.  Pres  term  and 
poststonn  sea  surface  height  (SSH,  m)  from  altimetry  and  models  are  superposed  on  these  plots  as 
contours.  TMI  SSTs  are  three -day  averages  ending  September  10  and  16,  respectively.  Contour  interval 
for  SST  is  1°C,  SSH  is  0.1  m  and  ASST  is  I°C.  The  storm  track  of  Hurricane  Ivan  compiled  by  the 
National  Hurricane  Center  is  plotted  for  every  6-hour  interval  The  solid  box  (2°  x  2°,  8S°-86°W,  24° - 
26°N)  represents  the  area  of  interest  for  the  analyses  presented  here.  Ivan  passed  over  this  region  on 
15  September  2004.  Altimeter  indicates  a  warm  core  eddy  (WCE)  near  89.2°  W,  26,7°N  and  a  smaller 
cold  core  eddy  (CCE)  near  87.I°W,  24.9°N,  southeast  of  the  larger  WCE. 


advection  by  wind -driven  currents  was  an  important  mech¬ 
anism  for  the  mixed-layer  heat  balance. 

[3i]  While  SST  cooling  varied  by  nearly  —Q.5°C  between 
AF  and  FF,  there  were  significant  differences  between  FA 
and  FF.  The  larger  simulated  SST  cooling  in  FA  were  not  a 
result  of  enhanced  shear-induced  vertical  mixing  because 
storm  induced  currents  in  all  cases  were  nearly  the  same 


(Figure  7).  The  different  response  of  the  SST  in  the  three 
experiments  (AF,  FF  and  FA)  can  be  explained  in  part  by 
differences  in  the  prestomi  thermal  stratification*  A  shal¬ 
lower  (deeper)  Lhermocline  is  likely  to  produce  stronger 
(weaker)  upper-ocean  cooling  with  a  given  surface  wind 
force.  The  depth  of  20°C  isotherm  (D2o.  measure  of 
the nnoc line  depth)  prior  to  the  passage  of  Ivan  indicated  a 
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Figure  7,  Comparison  of  ocean  surface  currents  (m  s“l)  before  (10  September,  top)  and  during 
(15  September,  bottom)  the  storm  from  (left)  AF,  (middle)  FF  and  (right)  FA,  respectively.  Shaded 
regions  indicate  the  corresponding  depth  of  20°C  isotherm  (D20>  m)  except  poststorm  (16  September, 
bottom)  D2<j  is  shown. 


shallow  thermoeline  west  of  ^B7°W  in  all  experiments 
(Figure  7).  The  prestorm  (September  10)  values  of  D20  in 
AF  and  FF  remained  similar  (~S0  m),  while  there  was  a 
significantly  shallower  D20  by  40  m  in  FA,  It  is  these 
differences  in  D20  that  lead  to  different  SST  responses 
during  the  passage  of  Ivan,  Thus,  the  wind-driven  mixing 
associated  with  the  storm  produced  a  greater  cooling  in  FA 
where  the  prestorm  thermodine  was  shallower.  The  phys¬ 
ical  processes  leading  to  different  SST  responses  in  these 
experiments  are  quantified  in  the  following  sections. 

4,  Mixed  Layer  Heat  Budget 

[32]  What  were  the  physical  processes  affecting  upper- 
ocean  cooling  during  Ivan's  passage?  To  address  this,  wc 


calculate  the  mixed  layer  heat-budget  terms  from  these  runs 
and  will  present  them  in  the  following  sections  and  sum¬ 
marize  them  in  Table  2.  The  mixed- layer  processes  being 
discussed  here  were  taken  at  a  constant  depth  of  20  m 
(slab  =  20  m)  rather  than  lime  varying  MLD.  The  heai- 
budget  terms  are  also  presented  below  the  mixed-layer  for 
two  depths  (50  and  100  m)  to  determine  the  underlying 
processes  leading  to  subsurface  temperature  changes.  These 
terms  are  averaged  over  the  area  between  88°W  and  86°W 
and  24°N  and  26°N  (2°  x  2°  box,  Figure  6)  where  the 
maximum  upper-ocean  cooling  occurred.  The  heat-budget 
terms  can  be  written  as  QT  =  -  Qu+V  —  Qw  +  Qs  +  Qdv  + 
Qdh;  that  is,  the  rate  of  change  of  heat  storage  (Or)  is  the 
sum  of  horizontal  (Qu+v)>  vertical  (Qw)  adveclton*  surface 
heat  flux  (Qs),  vertical  (Qt>v)  and  horizontal  diffusion 


Table  2 ,  Prestorm  (14  September),  Storm  (15  September)  and  Poststorm  (16  September)  Heat-Budget  Terms  From  Various  Model 
Experiments  and  Percentage  Contribution  of  Each  Term  in  Parentheses'1 


Experiments 

14  Sep 

15  Sep 

16  Sep 

AF 

FF 

FA 

AF 

FF 

FA 

FF-A9 

FA-F9 

AF 

FF 

FA 

Qr 

-325 

-152 

157 

-2055 

-1421 

-2870 

-2839 

—  1824 

-306 

305 

2061 

Qs 

-177 

—  1 68 

-65 

-124 

-34 

256 

366 

59 

11 

137 

334 

(54.4) 

(110.5) 

(-41.2) 

(6.1) 

(2.4) 

(-8,9) 

(-12.9) 

(-3.2) 

(-3.7) 

(44.8) 

(16.2) 

Qnv 

—82 

-95 

-165 

-896 

-1045 

-2002 

-2095 

-1079 

-150 

-250 

“259 

(25,0) 

(62.8) 

(-104.9) 

(43.6) 

(73.6) 

(69.8) 

(73,8) 

(59.2) 

(48.9) 

(-81.8) 

(-12.6) 

Qlnv 

-47 

94 

204 

-964 

-229 

-4C0 

-443 

-629 

-178 

235 

1846 

05.0) 

(-62,1) 

(129.9) 

(47.0) 

U6.t) 

(16,0) 

(15,6) 

(34.5) 

(58,0) 

(77.0) 

($9.6) 

Qw 

-20 

16 

181 

-11 

—  1 14 

-666 

-669 

—  176 

It 

183 

140 

(6.2) 

(-10.8) 

(115  4) 

(3,5) 

(8.0) 

(23.2) 

(23.6) 

(9,6) 

(-3.5) 

(59.8) 

(6,8) 

Tlcat-budgcl  terms  from  experiments  FF-A9  and  FA-F9  are  included  only  for  1 5  September,  Unit  is  W  m  These  are  daily  snapshots  averaged  for  Ihe 
2q  x  2*  box  described  in  the  text. 
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Figure  8*  Area-averaged  (8S°— 86°  W,  24°— 26°N)  heat-budget  terms  in  W  m  2  at  three  depths  20  m, 
50  m  and  100  m  from  FF  (left,  top  three  panels),  AF  (middle  panels),  FA  (right,  top  three  panels)  and 
heat-budget  terms  at  20  m  from  FF-F9  (left,  last  row)  and  FA-F9  (right,  last  row)  during  September  2004. 
The  heat-balance  terms  include  rate  of  change  of  heat  storage  (QT),  surface  heat  flux  (Qs),  vertical 
diffusion  (Qnv),  horizontal  diffusion  (Qr>H)  horizontal  advection  (Qu+V)  ailc*  vertical  advection  CQw); 
qt-  _  Qu+V  -  Qw  +  QlJV  +  Qdh  +  Qs,  Since  horizontal  diffusion  (Qdh)  is  negligibly  small,  it  is  not 
plotted  These  are  daily  snapshots.  Note  that  for  clarity  y-axis  has  different  divisional  lengths. 


(Qdh)  terms.  Since  horizontal  diffusion  (Qdh)  is  negligibly 
small  it  is  not  plotted.  The  time  scries  of  heat-budget 
terms  at  three  depths  20  m,  50  m  and  100  m  which  are 
representatives  of  the  mixed  layer,  mixed  layer  base  and 
upper-thennocline*  respectively,  from  FF,  AF,  and  FA 
are  displayed  in  Figure  8.  These  are  model  daily  snapshots, 

4.1,  Control  Run  FF  (Nona simulation) 

[rt]  Overall,  prestorm  SST  changes  were  caused  primar¬ 
ily  by  the  surface  beat  flux  term  (Qs)  while  horizontal 
advection  (Qu+v)  contributed  to  the  posts  term  net  warming. 
Vertical  diffusion  (QDV*  wind-driven  mixing)  dominated 
(—1045  W  m-2,  74%)  the  net  upper-ocean  cooling 
(-1421  W  m^2)  during  the  passage  of  Ivan  and  horizontal 
advection  caused  a  loss  of  228,5  W  rcT 2  (16%)  of  heat  from 
the  ocean.  Following  the  passage  of  Ivan  on  September  16, 
the  upper-ocean  started  wanning  at  a  rate  304  W  in'”2  due 
to  die  combined  effects  of  surface  heal  flux  and  advection 
terms  (Qu+v+\v)  which  overwhelmed  the  net  cooling 
(-250  W  m_fc)  due  to  vertical  mixing.  For  die  next 
three  days  the  ocean  gained  an  average  664  W  m  2  of  heat 


through  surface  heat  flux  (149  W  m-2)  horizontal  advection 
(294  W  m"2)  and  vertical  advection  (QWp  221  W  m^2). 
Thereafter,  the  upper-ocean  continued  to  warm  at  ^225  W  irr2. 
An  important  result  consistent  with  earlier  studies  is  dial  die 
vertical  mixing  dominated  the  mixed-layer  heat-balance 
during  the  storm. 

[34]  A  maj or  di  fferencc  in  heat  balance  terms  at  50  m  depth 
from  die  surface  (20  m)  was  the  role  played  by  vertical 
advection  in  wanning  or  cooling  at  the  mixed-layer  base. 
Throughout  the  period  the  rate  of  heat -storage  (Qx)  was  in 
general  caused  by  vertical  advection.  During  the  storm, 
upwelling  of  colder  water  from  below  (vertical  advection) 
was  responsible  for— 59%  of  cooling  (-2262  W  m  2)  while 
wind-driven  vertical  mixing  accounted  for  ^23%  of  cooling. 
Bodi  advective  terms  contributed  to  die  net  warming  after  the 
passage  of  Ivan.  Vertical  advection  (>90%)  controlled  the 
temperature  changes  at  1 00  m.  The  large  amplitude  temper¬ 
ature  changes  that  occurred  below  the  mixed  layer  (50  and 
100  m)  during  die  poststorm  period  were  a  result  of  the 
inertial  oscillations  with  a  period  near  28  hours,  generated  by 
the  onset  of  hurricane-force  winds  (not  shown).  Shear- 
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Figure  9.  Box  averaged  three-hourly  snapshots  of  model  net  surface  heat  flux  (QshfX  latent  heat  flux 
(Qlat)*  sensible  heat  flux  (Qsen),  long-wave  radiation  (Qlwr)  and  shortwave  radiation  (Qswr)  in  W  m-2 
from  AF. 


induced  entrainment  mixing  forced  by  near- inertial  motions 
after  the  passage  of  Gilbert  was  also  noted  by  Jacob  et  al 
[2000]. 

4,2.  Experiment  AF  (Assimilation  Initial  Condition) 

[35]  A  major  difference  between  FF  and  AF  was  the 
location  of  WCE  (Figures  6  and  7)  which  eventually  altered 
ihe  location  of  maximum  cooling  in  AF  Thus  a  comparison 
of  heat-balance  terms  between  AF  and  FF  delineate  the 
influence  of  the  WCE  on  SST  cooling.  As  is  evident  in 
Figure  8,  ihe  strong  ocean  currents  associated  wiih  the  WCE 
in  AF  clearly  contributed  to  the  mixed-layer  heat  budget. 
The  presto  no  heat-storage  rate  agreed  reasonably  well  with 
FF.  The  upper-ocean  began  to  cool  on  September  1 1  due  to 
the  combined  effects  of  surface  heat  loss  from  the  ocean 
(Qs)  and  horizontal  advection  (Qu+v).  About  54%  of  the 
total  heat  loss  on  September  14  resulted  from  the  air-sea 
heat  flux  and  vertical  mixing  accounted  for  25%.  As 
expecied,  this  high  surface  heat  loss  was  primarily  due  to 
the  high  latent  heat  flux  loss  (-473  W  m-2)  caused  by 
strong  winds  (Figure  9).  The  upper  ocean  cooled  at  a  rate 
—2055  W  in  2  during  the  storm;  vertical  diffusion  (turbu¬ 
lent  wind-driven  mixing)  contributed  44%  of  total  heat  loss 
and  horizontal  divergence  (cooling)  accounted  for  47%. 
Sensible  heat  flux  during  the  storm  period  showed  a  heat 
gain  by  the  ocean  (90  W  rrT2)  due  to  wanner  air-temper- 
ature.  Heat  loss  due  to  longwave  radiation  was  only  10- 
20  W  m“\  A  significant  part  of  the  warming  that  occurred 
following  the  passage  of  Ivan  resulted  from  net  surface  heat 
gain  by  the  ocean.  The  net-heat  gain  by  the  ocean  during 
this  period  was  primarily  due  to  a  reduction  in  latent  heal 
flux  by  weaker  winds  and  cool  SST's  (Figure  9).  The  solar 
radiation  (Qswr)  showed  a  marked  increase  from  87  W  m“2 
(15  September)  to  240  W  m  2  (17  September).  Following 
this  warming,  horizontal  advection  caused  mixed  layer 
cooling  during  the  period  20-22  September  2004,  While 
47%  of  upper-ocean  cooling  in  AF  wras  due  to  the  horizontal 
advection,  it  was  only  16%  in  FF  suggesting  that  WCE  had 
an  impact  on  the  SST  cooling. 

[36]  In  addition  to  both  vertical  mixing  and  horizontal 
divergence,  vertical  advection  of  cold  water  associated  with 
the  storm  provided  further  cooling  at  50  m.  Consistent  with 
FF,  at  100  m  depth,  the  rate  of  change  of  temperature  was 
primarily  due  to  vertical  advection  with  a  small  contribution 
from  horizontal  advection.  The  important  result  here  is  that 


the  horizontal  advection  is  as  large  as  the  vertical  diffusion 
term  in  the  overall  heat  budget 

4.3.  Experiment  FA  (Nonassimilation  Initial  Condition) 

[37]  An  important  shortcoming  of  the  above  experiments 
was  the  inconsistency  in  the  location  of  the  LC  and  WCE 
with  observations.  This  deficiency  was  overcome  in  part  by 
assimilating  SSH  on  to  the  model.  The  heat-budget  terms 
from  FA  were  significantly  different  from  the  two  other  frec- 
runs  especially  the  amplitude  of  the  wanning  and  cooling 
(Figure  8).  Except  during  Ivan's  passing,  all  ihe  terms 
showed  a  net  warming  tendency  in  die  overall  heat  budget. 
Vertical  mixing  dominated  the  mixed  layer  cooling  during 
Ivan’s  passage  that  accounted  for  70%  of  the  heat  loss 
which  was  comparable  with  that  in  FF  (74%),  The  contri¬ 
bution  from  surface  heat  flux  resulted  in  a  net  warming, 
primarily  due  to  the  large  reduction  of  latent  heat  flux 
(M8%  from  AF)  associated  with  excessive  SST  cooling. 
Unlike  the  other  two  free-runs,  in  FA,  there  was  an 
unexpected  increase  in  the  vertical -advection  contribution 
(23%)  to  the  total  cooling.  Since  these  experiments  arc 
forced  with  identical  wind  fields,  differences  in  the  up  well¬ 
ing  driven  by  the  curl  of  the  wind  stress  arc  very  unlikely. 
The  physical  processes  leading  to  these  differences  are 
further  discussed  in  the  next  section  5,  Horizontal  advection 
caused  16%  of  the  cooling  and  the  air-sea  exchange 
accountable  for  the  remainder  (9%  warming)  during  the 
storm,  A  comparison  of  balance  terms  between  FA  and  FF 
suggests  that  the  magnitude  of  the  cooling  actually  in¬ 
creased  dramatically  in  FA  because  of  an  increase  in  vertical 
advection  (or  upwelling)  during  the  storm.  Horizontal 
advection  dominated  the  heat  balance  during  16-17  Sep¬ 
tember  2004  following  Ivan’s  passage.  During  this  period, 
horizontal  advection  was  responsible  for  89%  of  upper- 
ocean  warming  and  ^14%  came  from  the  surface  heat  flux 
term.  Both  advection  terms  contributed  to  the  net  upper- 
ocean  wanning  for  several  days  after  18  September  2004, 
[.ik]  Beneath  die  mixed  layer  (at  50  and  100  m),  temper¬ 
ature  was  changed  primarily  by  vertical  advection.  Al  50  m 
depth  the  sum  of  vertical  (65%)  and  horizontal  (22%) 
advective  terms  caused  a  net  cooling  and  vertical  diffusion 
accounted  for  the  remainder.  Wind-driven  vertical  mixing 
was  still  a  significant  contributor  to  the  upper-ocean  cooling 
at  this  depth.  Upwelling  was  accountable  for  ^92%  of  the 
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Figure  ID*  (left)  Model  SST  (shaded,  °C)  and  SSH  (contours,  m)  on  1  September  from  FF  and  FA, 
(right)  Same  fields  from  FF-A9  and  FA-F9,  Contour  interval  tor  SST  (SSH)  is  1°C  (0.1  m). 


total  -15558  W  m  2  cooling  at  100  m  depth  during  the 
storm. 

4.4,  Experiments  FA-F9  and  FF-A9 
[w]  While  die  assimilation  of  SSH  in  FA  improved  the 
simulated  location  of  LC  and  WCE,  the  inability  to  get  the 
observed  amplitude  of  SST  cooling  during  the  storm 
remained  a  major  issue.  FA-F9  and  FF-A9  are  designed 
(Table  1)  to  unravel  this  issue  and  to  identify  the  scheme 
lhat  most  realistically  simulates  the  storm- induced  upper- 
ocean  response.  The  initial  (September  1)  and  final 
(September  16)  SST  (mixed  layer  temperature)  and  SSH 
from  these  runs  arc  shown  in  Figure  10.  The  initial  SST's 
(September  1 )  in  FF  and  FA  were  >29°C  almost  everywhere 
in  the  GoM  despite  the  marked  differences  in  the  location  of 
the  WCE.  The  cooling  in  FF-A9  due  to  Ivan’s  passage  was 
in  many  ways  identical  to  that  in  FA  (strong  cooling  in  the 
south  and  northeast  of  the  WCE)  though  cooling  of  the 
latter  area  was  less  pronounced  in  FA  (Figure  6).  The  area  of 
extreme  cooling  in  FA-F9  shifted  farther  south  (north)  with 
respect  to  FF  (AF),  analogous  with  the  observations.  The 
WCE  drifted  further  southweslward  from  its  initial  position 
in  the  free-n inning  model  whereas  the  location  of  the  WCE 
during  the  16  day  period  was  mostly  unchanged  in  the 
assimilative  run.  Two  conclusions  can  be  drawn  here  (1 )  the 
storm-induced  changes  in  the  mixed  layer  were  relatively 
insensitive  to  the  prestorm  conditions  when  assimilation 
was  included  and  (2)  tree-run  initialized  with  assimilation- 
model  fields  a  few  days  prior  to  the  storm  reproduced  the 
observed  changes  in  the  upper-ocean  reasonably  well, 

[40]  To  quantify  the  physical  processes  leading  to 
the  differences  in  the  mixed  layer  temperature  we  show 
hcat-baJanec  terms  at  20  m  depth  from  these  experiments 
(Figure  8),  A  comparison  of  heat  budget  terms  between  FF- 


A9  and  FA  delineate  the  sensitivity  of  the  model  to  initial 
conditions  when  assimilation  is  included.  The  adjustment  of 
the  model  initial-state  to  the  assimilative  fields  (SSH)  was 
indicated  by  a  sharp  jump  in  balance  terms  during  the  first 
7  days  in  FF-A9,  the  amplitude  of  the  warming  trend  in 
FF-A9  before  the  passage  of  I  van  was  much  smaller  than  that 
in  FA,  Regardless  of  these  differences,  the  physical  processes 
leading  to  maximum  cooling  during  Ivan's  passage  in  both 
experiments  (FA  and  FF-A9)  were  quite  similar.  The  vertical 
diJTusioncontribution  (mixed  layer  cooling)  that  increased  by 
4%  in  FF-A9  was  balanced  by  the  warming  caused  by  surface 
heat  flux  (Qs).  Due  to  the  differences  in  the  location  of  WCE 
the  rate  of  cooling  in  FA-F9  increased  by  403  W  m"2  from 
FF.  These  different  cooling  rates  were  caused  primarily  by 
differences  in  the  horizontal  advection  which  increased  from 
16%  (—229  W  itT2,  FF)  to  35%  (-629  W  nr3,  FA-F9).  The 
differences  in  the  vertical  mixing  between  the  two  runs 
accounted  for  34  W  uT2  cooling  in  FA-F9.  Horizontal 
advective  cooling  was  also  found  important  in  AF  (Figure  8) 
during  the  storm  period.  The  different  cooling  rates  between 
AF  (-2055  W  m~2)  and  FA-F9  (—1824  W  m~2)  were  caused 
in  part  by  differences  in  horizontal  advection  owing  to  differ¬ 
ences  in  the  WCE  location.  Thus,  depending  upon  the  location 
of  the  WCE  and  LC  region  horizontal  advection  contributed 
between  ^  1 6%  and  47%  to  the  total  mixed-layer  heal  budget. 

5*  Discussions 

[4i]  The  impact  of  data  assimilation  on  the  storm-induced 
response  during  die  passage  of  Hurricane  Ivan  (2-26  September 
2004)  is  examined  in  a  nested  regional  Gulf  of  Mexico 
(GoM)  HYCOM  model.  This  is  achieved  by  two 
approaches:  In  the  first  approach  (AF),  we  initialize  the 
model  with  assimilation  fields  that  include  loop  current 
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(LC)  and  warm  core  eddies  (WCE)  in  the  correct  locations. 
In  Hie  second  approach  (FA),  the  assimilation  is  introduced 
during  the  model  run  time.  Comparison  of  these  two  cases 
with  a  simulation  that  does  not  include  data  assimilation 
(FF)  delineates  the  effects  of  data  assimilation.  The  upper- 
ocean  responses  simulated  by  these  three  experiments  (AF, 
FA  and  FF)  were  similar,  but  there  were  large  differences  in 
the  location  and  amplitude  of  upper-ocean  cooling.  The 
differences  in  the  location  can  be  related  to  the  differences 
in  the  simulated  location  of  the  WCE  by  these  experiments. 
The  different  rale  of  SST  cooling  in  these  three  experiments 
can  be  explained  in  part  by  differences  in  the  prestorm 
thermal  stratification.  A  shallower  (deeper)  tbermocline 
together  with  a  stronger  (weaker)  vertical  temperature 
gradient  is  likely  to  produce  stronger  (weaker)  upper-ocean 
cooling  with  a  given  surface  wind  force.  The  depth  of  the 
2G°C  isotherm  (Djo,  measure  of  thermocline  depth)  prior  to 
the  passage  of  Ivan  indicated  a  shallow  thermocline  west  of 
~87°W  in  all  experiments  (Figure  7),  However,  D2o  in  FA 
(^40  m)  was  much  shallower  than  that  in  AH  and  FF 
(—80  m).  Also  noteworthy  is  the  larger  areal  extent  of  this 
shallow  zone  in  FA.  The  SSH  in  this  region  (west  of  LC) 
indicated  the  signature  of  a  cyclonic  eddy  which  appears  to 
be  generated  by  the  shear  instability  associated  with  the  LC 
(Figure  6).  While  a  shallower  D2o  is  expected  because  of 
this  cyclonic  eddy,  the  strength  of  the  eddy  simulated  by  AF 
and  FF  are  weaker  than  FA.  This  results  in  a  deeper  D2o  in 
AF  and  FF  in  comparison  with  FA. 

[42]  The  enhanced  SST  cooling  in  those  runs  that 
included  data  assimilation  were  not  a  result  of  differences 
in  vertical  mixing  or  upwelling  from  the  nonassimilative 
runs  because  under  identical  wind-forcing,  they  ought  to  be 
the  same.  A  possible  candidate  for  this  enhanced  SST 
cooling  is  the  difference  in  the  upper  thermocline  temper- 
amre  gradient.  To  explore  this  further,  we  show  in  Figure  1 1 
the  time  series  of  box  average  vertical  temperature  gradient 
in  the  upper  200  m  from  FF,  FA,  FF-A9  and  FA-F9, 
Included  also  in  Figure  1 1  is  ihc  mixed  layer  depth.  Both 
model  runs  that  included  data  assimilation  (FA  and  FF-A9) 
indicated  a  thin  prestorm  mixed  layer  and  a  sharp  upper 
thermocline  temperature  gradient  which  in  turn  enhanced 
ihc  upper-ocean  cooling  via  vertical  mixing  and  upwelling. 
On  the  other  hand,  a  weak  vertical  temperature  gradient  and 
somewhat  thicker  mixed  layer  during  the  prestorm  period  in 
FF  and  FA-F9  limited  the  cooling.  Vertical  sections  of 
temperature  (September  16),  vertical  diffusion  coefficient 
(KTi  m2  s  September  15)  and  vertical  velocity  (10  4  m  s 
September  15)  along  24°N  arc  depicted  in  Figure  11.  The 
location  of  maximum  vertical  mixing  occurred  to  the  right 
of  Ivan’s  path  as  expected  because  the  area  of  strongest 
wind  speeds  in  a  moving  hurricane  is  in  the  northeast 
quadrant.  The  deep  thermocline  (and  high  heat  storage) 
associated  with  the  LC  east  of  the  track  greatly  limited  the 
availability  of  cold  water  to  the  surface  from  below.  In 
contrast,  the  coldest  surface  water  was  located  somewhat  to 
the  left  of  Ivan’s  track. 

[43]  The  coincidence  of  the  maximum  cooling  with  the 
shallowest  thermocline  depth  (20°C  isotherm  deplh  was 
located  at  —75  m  or  less)  in  all  experiments  during  Ivan’s 
passage  suggests  that  both  upwelling  and  vertical  mixing 
contributed  to  the  upper-ocean  cooling.  It  should  be  noted 
that  the  shallow  thermocline  (88°-*86°W)  in  Figure  11 


(after  Ivan)  below  —100  m  were  not  greatly  different  from 
those  before  the  storm.  If  upwelling  occurs  simultaneously 
with  vertical  mixing,  then  the  mixed  layer  becomes  shal¬ 
lower  which  in  turn  enhances  surface  cooling.  The  vertical 
velocity  of  the  assimilative  runs  (FA  and  FF-A9)  predicted 
slightly  higher  values  (—15  x  10“4  m  s’1)  than  the  free- 
runs  (—10  x  10  4  m  s_l).  The  center  of  the  upwelling  was 
closer  to  the  center  of  the  storm  and  there  was  no  bias  to 
the  right  of  the  storm  track.  A  further  increase  in  surface 
cooling  is  possible  if  the  upper  thermocline  temperature 
gradient  is  large.  The  different  verlicaLmixing  cooling  rates 
between  FF  (-1045  W  m  2)  and  FA  (-2002  W  m  2)  were 
caused  primarily  by  differences  in  the  initial  thermal 
stratification  (Figure  1 1 )  and  could  not  be  explained  by 
the  differences  in  the  vertical  diffusion  coefficients.  For 
example,  the  volume  average  values  of  Kr  (880-86°W, 
24°-26°N,  0-50  m)  from  FF  and  FA  on  September  15 
was  0.044  and  0,028  nrf  s~\  respectively,  which  did  not 
support  a  greater  mixing  across  the  mixed  layer  base  in 
FA.  Thus,  it  is  evident  that  the  shallower  thermocline 
along  with  stronger  upwelling  and  a  larger  vertical 
temperature  gradient  in  the  assimilation  runs  (FA  and 
FF-A9)  resulted  in  much  stronger  upper-ocean  cooling 
than  the  runs  that  do  not  include  data  assimilation  (AF, 
FF  and  FA-F9). 

[44]  Allhough  the  model  results  presented  here  captured 
the  basic  details  of  the  processes  leading  to  die  upper-ocean 
cooling  response  to  Ivan,  there  are  several  model-data 
discrepancies  that  need  to  be  addressed.  The  experiments 
reported  here  (nonassimilative  runs)  dramatically  illustrate 
the  sensitivity  of  the  timing  of  eddy-shedding  behavior  to 
changes  in  the  initial  conditions.  This  issue  will  be  explored 
fully  in  the  near  future  through  a  scries  of  model  sensitivily 
experiments.  Though  surface  heat  fluxes  are  found  insig¬ 
nificant  in  the  mixed-layer  heat  budget,  considerable  uncer¬ 
tainty  remains  in  the  transfer  coefficients  that  need  further 
investigation.  The  results  presented  here  are  likely  sensitive 
to  the  choice  of  vertical  mixing  scheme  (C3SS)  employed 
here.  However,  the  overall  rate  of  cooling  and  mixed-layer 
deepening  would  remain  qualitatively  similar.  Finally,  a 
shallower  thermocline  and  an  associated  excessive  SST 
cooling  in  the  assimilative  runs  remain  a  concern.  This 
deficiency  of  the  model  may  be  associated  with  the  Cooper 
and  Haines's  [1996]  downward  projection  technique  that 
was  used  here;  this  issue  will  be  addressed  in  the  near 
future. 

6,  Summary 

[45]  A  20  layer  in  the  vertical  and  1/25°  (—4  km) 
horizontal  resolution  HYCOM  Gulf  of  Mexico  (GoM) 
model  was  employed  to  investigate  the  three-dimensional 
oceanic  processes  responding  to  Hurricane  Ivan  which 
traversed  the  GoM  during  the  period  14-16  September 
2004,  By  conducting  several  model  experiments,  we  exam¬ 
ined  die  sensitivity  of  the  location  of  the  WCE  and  die  LC 
to  the  initial  conditions  and  compared  a  nonassimilative 
simulation  with  one  that  included  assimilation  of  SSH.  The 
storm  induced  changes  were  insensitive  to  the  initial  con¬ 
ditions  in  the  assimilative  runs  but  highly  sensitive  in  the 
nonassimilation  owing  to  variations  in  the  trajectories  of  the 
WCE.  SST  decreased  — 6°C  as  Ivan  passed  over  the  GoM  a! 
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Figure  11,  (top)  Time  scries  of  vertical  temperature  gradient  (°C  m”\  contours)  and  mixed-layer  depth 
(m,  doited  line)  averaged  over  88D-860W,  24°-26°N.  (middle)  Meridional  vertical  sections  (24°N)  of 
temperature  (contours)  on  16  September  (dashed  line  indicates  the  climatological  September  mean 
position  of  the  depth  of  20°C  isotherm  from  the  General  Digital  Environmental  Model  (GDEM),  and 
shaded  regions  indicate  vertical  diffusion  coefficient  (ln(Xr]t  natural  logarithm  of  Kr  in  m2  s-1)), 
(bottom)  Vertical  velocity  (>5  x  IG~4  m  stI)  from  FF,  FA,  FF-A9  and  FA-F9  on  15  September  2004. 
Longitudinal  position  of  Ivan's  eye  at  24 °N  is  marked. 


6.5  m  s_l  in  agreement  with  the  TMI  observations.  The 
location  of  maximum  cooling  associated  with  Ivan's  pas¬ 
sage  varied  depending  upon  the  location  of  the  WOE. 
However,  maximum  cooling  in  all  experiments  occurred 
at  the  center  of  the  storm  track  (south  of  ^26°N);  this 
differs  from  several  observations  and  modeling  studies 
where  cold  water  areas  were  found  to  the  right  of  the  storm 
track.  While  the  maximum  vertical  mixing  was  experienced 


in  the  northeast  quadrant  of  the  storm,  the  lack  of 
corresponding  SST  response  was  due  to  the  underlying 
thermal  structure  of  the  water  column.  Inclusion  of  data 
assimilation  did  not  greatly  alter  the  pattern  of  SST  cooling, 
but  caused  pronounced  SST  cooling  (>8°C)  over  a  larger 
area.  In  the  WCE  region  the  SST  cooling  wras  only  1°-2°C 
[46]  The  physical  processes  underlying  the  uppcr-occan 
cooling  that  occurred  during  the  hurricane  passage  were 
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quantified.  Overall,  -64%  of  the  cooling  was  due  to 
vertical  mixing  caused  by  the  turbulence  generated  from 
the  strong  shear-stress  across  the  base  of  the  mixed  layer. 
However,  this  value  varied  from  58.8%  (all  nonassimilative 
runs)  to  71.8%  (all  assimilative  runs)  depending  on  the 
model  experiments.  In  all  experiments,  die  vertical  d illu¬ 
sion  term  contributed  to  the  mixed-layer  heat  budget  only 
during  13  -17  September.  The  dominant  role  of  entrainment 
mixing  was  also  previously  noted  by  Jacob  et  al  [2000] 
based  on  observations  and  by  Price  [1981]  based  on  model 
results.  Vertical  advection  (upwelling)  caused  a  significant 
portion  of  cooling  (23.4%)  in  those  runs  that  included  data 
assimilation;  a  three  fold  increase  from  the  free  runs 
average  (7%).  This  enhanced  upper-ocean  cooling  was 
caused  primarily  by  the  initial  thermal  stratification;  a 
shallower  thermocline  (Figure  7)  and  a  stronger  upper- 
thermodiue  gradient  (Figure  11).  In  all  experiments  the 
air-sea  exchange  was  a  small  component  of  the  mixed  layer 
heat  budget  which  overall  accounted  for  —4%  and  in  some 
cases  actually  contributed  to  net  surface  warming.  The 
mixed-layer  heat  balance  following  the  passage  of  Ivan 
was  dominated  by  horizontal  advection  due  to  oscillations 
in  velocity  components  which  are  associated  with  inertia- 
gravity  waves  induced  by  the  storm.  Vertical  velocities  with 
alternating  upward  and  downward  oscillations  associated 
with  the  inertial-gravity  waves  also  persisted  in  the  wake  of 
the  storm.  Overall  the  rale  of  cooling  below  the  mixed  layer 
(50  m)  was  largely  due  to  vertical  divergence  (—54%);  the 
combined  effects  of  horizontal  advection  (—26%)  and 
vertical  mixing  (—20%)  accounted  for  the  remainder.  The 
storm  induced  vertical  motion  led  to  large  vertical  advec¬ 
tion  at  100  m  which  was  responsible  for  —90%  of  the  rate 
of  heat  change, 
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